Quantum ripples in strongly correlated metals 
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We study how well-known effects of the long-ranged Friedel oscillations are affected by strong 
electronic correlations. We first show that their range and amplitude are significantly suppressed 
in strongly renormalized Fermi liquids. We then investigate the interplay of elastic and inelastic 
scattering in the presence of these oscillations. In the singular case of two-dimensional systems, we 
show how the anomalous ballistic scattering rate is confined to a very restricted temperature range 
even for moderate correlations. In general, our analytical results indicate that a prominent role of 
Friedel oscillations is relegated to weakly interacting systems. 
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Introduction. — Fermi liquid theory is known to success- 
fully describe the leading low temperature behavior of 
metals, even in instances of very strong correlations (e.g. 
heavy fermions [lfl). In the presence of perturbations that 
break translational symmetry, such as impurities and de- 
fects, the Fermi liquid re-adjusts itself producing a spa- 
tially inhomogeneous pseudo-potential "seen" by quasi- 
particles 0, [J Here the wave nature of the electrons is 
manifested by the formation of "ripples", the Friedel os- 
cillations [4|,l5j], surrounding the perturbation. Scattering 
processes of quasiparticles off these ripples then produce 
new corrections to the T-dependence in transport quan- 
tities Q. 

How significant are these corrections? The answer, 
of course, depends on how broad the dynamic range is, 
in which such leading order non-analyticities dominate. 
This question, as usual, cannot be answered by Fermi liq- 
uid theory itself. What is needed is a microscopic model 
calculation that is not restricted to obtaining the form of 
leading terms. A careful and precise model calculation 
with such a goal is the central content of this paper. 

We focus on single non-magnetic impurity scattering in 
an otherwise uniform strongly interacting paramagnetic 
metal, where the analysis is most straightforward and 
transparent, but this general issue is of key relevance also 
for the diffusive regime. Our mostly analytical results 
demonstrate that: (i) for sufficiently weak correlations 
we recover the results of the Hartree-Fock approximation, 
in which the effective scattering potential generated by 
the impurity is set by the long-ranged Friedel oscillations; 
(ii) as we approach the Mott transition, however, these 
oscillations are strongly suppressed as the charge screen- 
ing becomes more and more local, corresponding to a 
shorter "healing length"; (iii) a combination of "healing" 
and inelastic scattering strongly suppresses the Friedel 
oscillation effects even for moderate correlations. 

Model and method. — We study the paramagnetic phase 
of the disordered Hubbard model on a cubic lattice in d 
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Figure 1: (a) Electronic density deviations 5ru — m — 1 (see 
text) displaying characteristic Friedel oscillations. From top 
to bottom we have m/m* = 1.00 and 0.30. The Friedel os- 
cillations appear here as crosses because of the underlying 
Fermi surface anisotropy [?[. As we enter the strongly corre- 
lated regime, these oscillations are suppressed, similar to the 
"healing" effect found in Ref. [sj. (b) Quasiparticle weight 
deviations SZi = Zi — Zq. From top to bottom we have 
m/m* — 0.85 and 0.08. While for moderate values of inter- 
actions SZi also displays Friedel oscillations, it has a leading 
exponential decay close to the Mott transition. Here, we have 
used a 30X30 square lattice with periodic boundary condi- 
tions and £ = —D. The color scales encode only the positive 
values of both Sm and SZi. 
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where tij are the hopping matrix elements between 
nearest- neighbor sites, c\ a (ci a ) is the creation (annihi- 
lation) operator of an electron with spin projection a at 
site i, U is the on-site Hubbard repulsion, — c\ a Ci a 
is the number operator, and £j are the site energies. 
All energies will be expressed in units of the clean half- 
bandwidth (Fermi energy) D and we approach the Mott 
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metal-insulator transition (MIT) by increasing U at half 
filling (chemical potential fi = U/2). To treat this 
model, we employ the slave boson mean-field theory of 
Kotliar and Ruckenstein which is equivalent to the 
Gutzwiller variational approximation [lOj. In this ap- 
proach, the renormalized site energies Vi and the lo- 
cal quasiparticle weights Zi are variationally calculated 
through the saddle-point solution of the cor resp onding 
Kotliar-Ruckenstein slave boson functional [llj. This 
theory is mathematically equivalent to a generalization 
of the dynamical mean field theory (DMFT) [12j to fi- 
nite dimensions, the statistical DMFT (statDMFT) 
implemented using a slave boson impurity solver |l4|. 
At T = and in the uniform limit (si = 0) we have 



Vq = and Zq = 1 — u l , with u = U/U c [^. The critical 
interaction value U c for which the Mott metal-insulator 
transition occurs is characterized by the divergence of the 
effective mass m* = m/Zo, where m is the electron band 
mass, indicating the transmutation of all electrons into 
localized magnetic moments. 

We consider first a generic weak disorder potential 
(|£j| <D) and expand the resulting mean- field equa- 
tions, Eqs. (5) and (6) from 14|, around the uniform 
solution. For particle-hole symmetry Zi = Zq + (ef) 
and, up to first order in £,-, the renormalized disorder 
potential, which is the effective potential "seen" by the 
quasiparticles at the Fermi level, reads (summation over 
repeated indices implied throughout) 



(2) 



The matrix M (u) is the lattice Fourier transform of 



M q (u) = 1 - 2g (u) 



(3) 



where n q ^ is the usual static Lindhard polarization func- 
tion fl5| . of the clean, non- interacting system and 

g (u) = (l + u)(l-u) 2 [2u + u 2 (1 - u)]^ 1 . (4) 

Charges rearrange themselves to screen the impurity 
potential and the local electronic density is given by rii = 
1 + Srii , where 

5m = -4g (u) [U c (1 - u 2 )] - 1 [M- 1 (u)].. e h (5) 

with the same spatial structure as Vi in Eq. @. 

We particularize now to a single impurity with energy 
s placed at the site o such that e,; = 5i^ e . Although 
obtained for weak scattering, our analytical theory does 
capture the qualitative trends even when the scattering 
is not weak, as we show numerically, see Fig. (fTJ). 

Weak and strong coupling limits. — In the weak cou- 
pling regime (u — > 0) there is no mass renormalization 
and Eqs. ([2]) and ((5|) agree with the Hartree-Fock so- 
lution of the Hubbard model, with a local (static) self- 
energy given by = Un®, where n° = 1 



the non-interacting electronic density and EEy is the lat- 
tice Fourier transform of TLq^ . Even though the bare 
impurity potential Si is localized in space, the density 
deviation 8n® displays long-ranged Friedel oscillations en- 
coded in Ely . For example, in a free electron gas we have 
5n (r) ~ cos (2kpr) /r d , where r is the distance to the im- 
purity and kp is the Fermi momentum. Slowly decaying 
Friedel oscillations are a direct consequence of the gap- 
less nature of Fermi liquid excitations. The renormalized 
disorder potential reads Vi ~ Si + UTl^Ej, implying that 
the electrons scatter not only off the local bare impurity, 
but also off the long-ranged potential generated by the 
Friedel oscillations. 

As we approach the critical region (u — > 1), however, 
the density deviation in Eq. ([5]) becomes 



Sn, = - 



2(1 -u) 

Uc 



— {1-uf 



n (°)l 



,(6) 



showing a suppression of the Friedel oscillations: the non- 
local part of 5rii is a factor (1 — it) 2 smaller than the local 
one. Therefore, the electronic density is significantly dis- 
turbed only in the vicinity of the impurities, implying 
a much shorter "healing length", see Fig. ([J) (a). The 
suppression of the slow spatial decay in Snt reflects the 
fundamental tendency of quasiparticles to become local- 
ized as the system approaches the Mott insulator. 
The renormalized disorder potential, Eq. @, is 



(1 



1 u c 1 [n<°>]^ 



(7) 



2ng ) e, is 



and the screened impurity potential is just as non-local 
as for small u, except for a reduction of the overall ampli- 
tude scale. Therefore, we should not be guided by den- 
sity fluctuations alone which are indeed healed very ef- 
fectively. However, we notice that Vi goes to zero linearly 
at all lattice sites at the transition, signaling a complete 
suppression of disorder by interactions (l6| . 

To obtain the leading energy correction of Zi, we have 
to expand the mean-field equations up to second order 
in £j. At intermediate values of the interaction, devia- 
tions in the quasiparticle weights SZ t = Zi — Z also show 
Friedel-like oscillations. Close to the Mott transition, 5Zi 
displays a leading exponential decay from the impurity 
site, since all additional terms describing its long-range 
oscillations are of higher order in (1 — u), see Fig. (TTJ) (b) 
and Eq. (J5|) below (the details of the calculation will be 
presented elsewhere). A finite impurity potential tends 
to push the site occupation away from half filling, thus 
reducing the tendency to form a local moment and ren- 
dering the given site locally more metallic by increasing 
Zi. As spatial correlations grow, this "metallization" of 
the correlated metal tends to spread out away from the 
impurity, thus creating metallic "puddles" (red region in 
Fig. (H}(b)) m an almost-localized host. A somewhat 
similar result emerges in the t—t'—J model, in which 
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an impurity induces a local staggered magnetization in 
its vicinity [l?], 18 1 whose spatial extent also increases 



with correlations. The critical behavior of SZi is captured 
by our analytical expressions and we can show that, for 
d > 2, and r io /£ > 1 



SZi 



1 -u 



_(l-d)/2 



JJ2 y 2 a+d)/2^(d-3)/2 



where r.,, 
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£ = (2z(l — u)) 1 plays the role of a correlation length. 
This correlation length diverges at the transition with a 
mean field exponent 1/2. Previous studies on the inter- 
face of a strongly correlated metal and a Mott insulator 
Til 20| . which use techniques similar to ours, also find an 
analogous leading exponential decay of the quasiparticle 
weight upon entering the Mott insulator from the metal. 
In those studies, however, the oscillating terms of Eq. f8"j) 
seem to have been overlooked. 

Leading finite T corrections for transport and inelas- 
tic cutoffs. — We would like now to go beyond T = and 
study the behavior of the leading temperature corrections 
to the resistivity p (T) as a function of the correlations. 
We focus henceforth only on 2d systems, since, in the 
weakly correlated regime, electron scattering by Friedel 
oscillations leads to a non-Fermi-liquid linear tempera- 
ture correction to p (T) [(| [2l[ in the ballistic regime. 

The transport scattering rate is given by 



rn<°>i 



(8) 



z is the coordination number, and 
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30) , (9) 



where Ui mp is the impurity concentration, T q is the 
Fourier transform of the T matrix, q — 2ksm0 /2 is the 
transferred momentum, 9 is the scattering angle. To sim- 
plify our analytical expressions, we henceforth focus on 
the case a free electron dispersion e = k 2 /2m; we care- 
fully checked that no significant changes are found if a 
different dispersion is used. Up to first order in the im- 
purity potential e D , the T matrix is simply given by the 
renormalized disorder potential (Eq. (J2J). The scatter- 
ing time is given by the average r^ r = f der^ x (e) /' (e), 
where /' (e) is the derivative of the Fermi distribution 
function. 

In our slave boson mean-field theory, the electronic 
self-energy is purely real [j| and describes only the elas- 
tic scattering of the electrons off a temperature depen- 
dent screened impurity potential. However, this scheme 
should really be regarded as a variational calculation of 
the quasiparticle parameters within our statDMFT pro- 
cedure. In a fuller treatment, there is also an imaginary 
part in the self-energy, reflecting inelastic effects. For the 
purposes of examining the leading perturbative effects of 
impurity scattering, the imaginary part can be computed 
in the uniform system, where it emerges naturally in the 



context of local Fermi liquid theories like DMFT [12j,[22| 
and is given by 



7 (T) = A (u) Tp (T/TpY 



(10) 



where Tp is the Fermi temperature and the function A (u) 
has the following limits: A (u) ~ u 2 for small u and 
A (it) ~ (1 — u)~ 2 ~ (m*/m) 2 close to the Mott transi- 
tion. These limits can be understood from the fact that 
in the weakly correlated regime inelastic scattering ef- 
fects are perturbative, whereas in the strongly correlated 
regime we recover the well known Kadowaki- Woods rela- 
tion , observed in several strongly correlated systems, 
and which holds within the DMFT picture we use. 

There are two leading contributions from inelastic scat- 
tering. A bulk one, present even in the clean limit, given 
^ T in (^) = VI (T), where r/ is a geometrical factor de- 
pending on the band structure used in the DMFT calcu- 
lation. t^ 1 (T) simply adds to r^ 1 (T) in Eq. © through 
Matthiessen's rule, since we consider very dilute impuri- 
ties. In addition to this, a finite imaginary part also cuts 
off the leading non-analyticities of the elastic scattering 
off Friedel oscillations 

[alii- This is taken into account 

in the calculation of riq -* , considering that the electron 
energy has now an imaginary part given by 7 (T) . The 
calculation of riq ' in the presence of inelastic broadening 
is carefully discussed in Ref. [23( for the 2d electron gas 
and we use the analytical form of n q as obtained there. 
The final form of r^ 1 (T), valid for T <C T>, reads 

r t - x (T) = r - 1 A 2 ( M )|l + 2^-«( U ) U ;(r, 7 (T))| 

+ rn(T), (11) 



where 
to (T, 7) 



T»* 3 (f) 



Re 



inr 



2tt 



7(T) 
T 



-ln(27r) 



1(T)_JT F 



2ttT V 2ttT 



(12) 



A(u) = g(u)[(l-u 2 )(p(e F )U c + g(u))] 1 , (13) 



a(u) = 2p(e F )U c [p(e F )U c + g{u)] 



(14) 



Tq 1 is the zero-temperature impurity scattering rate, 
p(sf) is the clean electronic density of states at the 
Fermi level and F (z) is the Gamma function. The func- 
tion A (u) controls the amplitude of the scattering rate 
from the screened impurity potential. In the weakly in- 
teracting regime A (u) ~ 1, whereas in the critical region 
A (u) ~ (1 — u) due to a vanishing V{ at U — U c . Eq. |7}. 
The function a (it) gives the strength of the leading tem- 
perature correction. It goes as U for weak correlations, 
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Figure 2: Limiting temperatures to the linear in T regime 
of the resistivity of a 2d electron liquid on a log-log scale, 
calculated under two different assumptions (see text) for two 
values of the parameter p(£f)U c [2J|. T* mp is almost inde- 
pendent of p(ef)U c and only one curve is shown here. T£ ulk 
is the dominating cutoff above which the linear behavior is 
lost. Here, we have used tqEf = 10, corresponding to a weak 
impurity potential. 



indicating that the temperature corrections only arise in 
the presence of electron-electron repulsion, and saturates 
to 2 close to the Mott transition. The function w (T, 7) 
encodes the dependence of the leading temperature cor- 
rection on 7 (T). 

For 7 (T) = 0, only elastic scattering is present and we 
obtain w (T, 0) = 0.5. Plugging this into Eq. (fTTj) , we see 
that the linear in T correction is found for all u > p| 
and is limited only by the overall amplitude A (u), which 
is in accordance with Eq. ([7]). 

However, for a finite 7(T), the linear region of p(T) 
considerably narrows as we enter the correlated regime. 
As there are two leading contributions from inelastic scat- 
tering, we analyze their individual effects separately by 
denning two threshold temperatures bounding the non- 
Fermi-liquid region from above. They are obtained by 
comparing the purely elastic t^, (T) (7 (T) — ► in Eq. 
(fTTj) ) with two other scattering rates: one in which we let 
77 — > {T* np ) and the other in which we set w (T, 7) — > 
w(T,0) (T* ulk ) in Eq. CP- We interpolate the func- 
tion A (u) using DMFT with both quantum Monte Carlo 
and iterated perturbation theory [25] as impurity solvers. 
We use rj = 10, a value also in agreement with DMFT 
calculations [25|. As we can see from Fig. ©, T£ ulk is 
strictly smaller than T* mp . This is because T£ ulk is not 
only proportional to the infinitesimal impurity concentra- 
tion riimp (or, equivalently, to (to£f) ) but, in the crit- 
ical region, T£ ulk ~ (m/m*f, whereas T* mp ~ {m/m*f . 
Thus, for any degree of correlations, it is T£ ulk which 
sets the upper bound on the linear in T region of p (T) . 
From Fig. ([2|) we see that, even for very moderate cor- 
relations, e.g. for m/m* ~ 0.9, T£ ulk ~ 10 -2 Tp, and for 
m/m* ~ 0.6, already T* ulk ~ 1Q- 4 T F . Thus, the non- 



Fermi liquid region is limited to very low temperatures. 
Ultimately, as the linear in T regime is also bounded from 
below by a crossover to the diffusive regime, the ballis- 
tic T-interval in which these elastic corrections dominate 
may not be present at all. 

Conclusions. — We presented a detailed, mostly ana- 
lytical, model calculation of the effects of a single non- 
magnetic impurity placed in a correlated host. We find 
that strong correlations tend to reduce the effects of the 
long-range part of the Friedel oscillations and our work 
provides clear analytical insight into how this happen. It 
should be possible to directly test our quantitative pre- 
dictions by means of current generation scanning tun- 
neling microscopy methods, shedding new light on the 
behavior near the Mott MIT. It is noteworthy that im- 
purities placed in d-wave superconductors also produce 
slowly decaying perturbations in real space, reflecting 
their gapless nature, through a mechanism closely re- 
lated to Friedel oscillations in normal metals. Recent 
work by Garg et al. Q shows, by using a very similar 
theoretical approach as we do, that in this system strong 
correlations also lead to spatial "healing". We believe 
that both phenomena have a closely related origin and 
our results strongly suggest that the "healing" effect is a 
more general property of correlated metals close to the 
Mott transition, not an effect specific to cuprates or the 
superconducting state. 
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